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•  An  analytical  method  for  optimal  condition  monitoring  thresholds  is  proposed. 

•  The  approach  regards  the  electrical  impedance  as  complex  random  variable. 

•  Impedance  statistical  properties  are  calculated  via  Morlet  wavelet  coefficients. 

•  Thresholds  are  determined  based  on  the  probability  of  false  alarm. 

•  Unified  and  scale-free  condition  indicator  for  PEM  fuel  cells  is  proposed. 
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When  commissioning  or  restarting  a  system  after  a  maintenance  action  there  is  a  need  to  properly  tune 
the  decision  thresholds  of  the  diagnostic  system.  Too  low  or  too  high  thresholds  may  implicate  either 
missed  alarms  or  false  alarm  rates.  This  paper  suggests  an  efficient  data-driven  approach  to  optimal 
setting  of  decision  thresholds  for  a  PEM  fuel  cell  system  based  solely  on  data  acquired  from  the  system  in 
reference  state  of  health  (i.e.  under  fault  free  operation).  The  only  design  parameter  is  the  desired  false 
alarm  rate.  Technically,  the  problem  reduces  to  analytically  determining  the  probability  distribution  of 
the  fuel  cell's  complex  impedance  and  its  particular  components.  Employing  pseudo-random  binary 
sequence  perturbation  signals,  the  distribution  of  the  impedance  is  estimated  through  the  complex 
wavelet  coefficients  of  the  fuel  cell  voltage  and  current.  The  approach  is  validated  on  a  PEM  fuel  cell 
system  subjected  to  various  faults. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Like  any  other  item  of  equipment  the  reliability  and  durability  of 
proton  exchange  membrane  (PEM)  fuel  cells  is  affected  by  faults  on 
various  components  [1,2].  These  faults  include  corrosion  of  the 
electrodes  and  degradation  of  membranes  [3—6],  catalyst  and 
membrane  poisoning  [7—10],  and  water  management  faults  (i.e. 
flooding  of  gas  channels  and  membrane  drying)  [11—13],  Timely 
detection  and  diagnosis  of  these  faults  is  of  great  importance  for 
future  employment  of  PEM  fuel  cell  technology  on  larger  scale 
[14.15], 

Majority  of  faults  that  may  occur  on  PEM  fuel  cells  can  be 
detected  with  multitude  of  approaches  based  on  electrochemical 
impedance  spectroscopy  (EIS)  [16—21],  EIS  characterization  of  fuel 
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cell  flooding  and  membrane  drying  was  reported  by  plethora  of 
authors  [22—28].  Portion  of  these  authors  used  EIS  measurement 
data  to  directly  detect  faults  [26—28],  meanwhile  others  proposed 
model  based  approaches  [22—25],  Employing  somehow  more 
advanced  approach,  Kadyk  et  al.  [29,30]  proposed  a  non-linear 
extension  of  the  EIS  diagnostic  approach  with  non-linear  fre¬ 
quency  response  analysis  (NFRA)  to  diagnose  fuel  cell  flooding  and 
membrane  dehydration.  In  addition  to  the  detection  of  water 
management  faults,  Le  Canut  et  al.  [27]  and  Kadyk  et  al.  [31]  also 
tackled  detection  of  anode  catalyst  poisoning  with  CO  by  employ¬ 
ing  more  complex  NFRA.  All  these  numerous  studies  confirm  that 
various  PEM  fuel  cell  faults  can  be  observed  through  EIS  impedance 
measurements.  Therefore,  an  autonomous  and  effectual  condition 
monitoring  system  for  fuel  cells  can  be  implemented  by  using  EIS 
features. 

To  the  best  of  the  authors'  knowledge,  at  this  stage  of  PEM  fuel 
cell  diagnostic's  development,  using  EIS  features  for  condition 
monitoring  requires  beforehand  EIS  characterisation  of  fuel  cell 
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behaviour  under  normal  and  faulty  operation.  Having  such  data, 
one  can  set  up  threshold  values  that  mark  the  border  between  fault 
free  and  faulty  regions. 

However,  performing  such  characterisation  procedures  may  be 
infeasible  in  many  practical  situations,  since  they  may  lead  to 
irreversible  damage  of  the  fuel  cell  system  (e.g.  extensive  drying 
and  CO  poisoning  [1,2]).  Moreover,  due  to  variation  in  parameters 
from  one  fuel  cell  to  another,  once  the  E1S  data  is  obtained  for  one 
particular  fuel  cell,  there  is  no  straightforward  way  to  transfer  this 
particular  knowledge  to  another  cell  without  conducting  new 
complete  set  of  characterisation  procedure.  Therefore  a  vital 
question  arises  how  to  design  a  condition  monitoring  system 
employing  measurements  conducted  solely  under  fault  free 
(normal)  operation  of  a  fuel  cell. 

Implementing  an  effective  condition  monitoring  system  can 
significantly  improve  the  operation  of  in-the-field  fuel  cell  power 
units  [32].  Many  of  the  problems  for  designing  an  effective  condi¬ 
tion  monitoring  system  have  been  addressed  in  the  context  of 
mechanical  systems  [33,34],  The  main  difficulty  is  to  determine  the 
statistical  properties  of  the  selected  features  and  model  their 
behaviour  under  faulty  conditions  [35,36],  Achieving  this  goal  has 
significant  practical  importance  for  any  system  integrator. 

Addressing  the  issue  of  beforehand  characterisation,  this  paper 
proposes  a  systematic  approach  for  determining  the  bounds  of  the 
fault  free  region  of  operation  without  inducing  any  faults.  In  the 
proposed  approach  the  fuel  cell  impedance  is  considered  as  a 
complex  random  variable.  The  corresponding  probability  distribu¬ 
tion  functions  are  derived,  allowing  complete  statistical  charac¬ 
terisation  of  the  behaviour  of  the  complex  impedance.  As  a  result, 
the  only  design  parameter  is  the  desired  false  alarm  rate,  which  is 
intuitive  and  veiy  well  understood  even  by  non-specialists  in  the 
area  of  condition  monitoring.  The  effectiveness  of  the  proposed 
approach  was  evaluated  on  a  8.5  kW  PEM  fuel  cell  subjected  to 
various  water  management  faults. 

The  first  step  towards  specifying  the  fault  free  region  is  the 
analysis  of  the  statistical  properties  of  the  feature  set  as  described 
in  Section  2.  Using  these  statistical  properties  the  method  for 
calculating  the  bounds  of  the  fault  free  region  as  well  as  the  overall 
condition  indicator  are  presented  in  Section  3.  Experimental  eval¬ 
uation  of  the  proposed  condition  indicator  is  presented  in  Section  4. 

2.  Statistical  properties  of  the  fuel  cell  impedance 

Fuel  cell  electrical  impedance  is  a  complex  variable  that  is 
defined  as  a  frequency  domain  ratio  of  the  voltage  to  the  current. 
Regarding  the  fuel  cell  system  as  linear  and  time  invariant  one,  the 
electrical  impedance  can  be  referred  as  its  transfer  function  which 
has  a  deterministic  nature  for  each  frequency.  However,  when 
dealing  with  real  world  scenarios  the  observed  fuel  cell  system  is 
influenced  by  external  disturbances,  which  are  generally  modelled 
as  stochastic  signals  [37],  as  shown  in  Fig.  1.  Since  these  distur¬ 
bances  are  in  many  cases  unmeasurable,  their  influence  is  incor¬ 
porated  into  the  model  thus  making  the  mapping  between  the 
current  and  voltage  a  stochastic  one.  Therefore,  the  observed 
electrical  impedance  should  be  treated  as  a  complex  random 
variable. 

disturbances 
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Fig.  1.  Representation  of  fuel  cell  system  excited  by  current  and  affected  by 
disturbances. 


In  order  to  properly  describe  the  electrical  impedance  as  a 
complex  random  variable,  the  first  step  is  determining  its  proba¬ 
bility  distribution  function.  For  that  purpose  the  fuel  cell  system 
can  be  excited  using  (PRBS)  signal  [38],  Following  the  procedure 
described  by  Debenjak  et  al.  [38  ,  the  impedance  characteristic  can 
be  efficiently  estimated  by  using  the  complex  wavelet  coefficients 
of  the  continuous  wavelet  transform  of  the  fuel  cell's  current  i(t) 
and  voltage  u(t)  using  the  Morlet  mother  wavelet,  which  is  used  in 
a  variety  of  fields  [38—41  ]. 

2.3.  Continuous  wavelet  transform  using  the  Morlet  wavelet 

The  continuous  wavelet  transform  (CWT)  of  a  square  integrable 
function /( t)  eL2(R)  is  defined  as  [42]: 

CO 

Wf(S,u)  =  (f(t),l/rus(t))  =  J  f(t)^s(t)dt,  (1) 

—  co 

where  ^u>s(t)  is  a  scaled  and  translated  version  of  the  mother 
wavelet  ^(t): 

(2) 

The  wavelet  coefficients  (1)  describe  the  analysed  signaler)  on 
the  time-scale  plane.  The  conversion  between  scale  s  and  fre¬ 
quency  /  is  straightforward  and  depends  on  the  selection  of  the 
mother  wavelet  (2).  Therefore,  in  the  remaining  of  the  paper 
instead  of  scale  s  all  subsequent  relations  will  rely  on  actual  fre¬ 
quency  /. 

For  the  impedance  estimation,  one  needs  information  about  the 
instantaneous  amplitude  and  phase  of  the  electrical  current  i(t )  and 
voltage  u(t).  Therefore,  the  Morlet  wavelet  was  chosen  as  a  mother 
wavelet  (2),  which  reads  [43]: 

y<t)  =  7r^a/4 -  e~wi>/2je~t2/2,  (3) 

where  ojq  is  the  ratio  between  the  highest  and  the  second  most 
highest  peak  and  is  usually  set  to  ojq  =  5. 

The  Morlet  wavelet  (3)  is  an  analytical  function  i.e.  it  has  only 
positive  frequencies.  As  a  result,  the  wavelet  coefficients  Wf[tf)  in 
( 1 )  are  complex  values  and  at  each  time  translation  t  and  frequency 
/  the  wavelet  coefficients  give  the  instantaneous  amplitude  and 
phase. 

2.2.  PRBS  as  excitation  signal 

PRBS  can  be  regarded  as  sufficiently  close  to  stationary 
random  noise  whose  amplitude  probability  distribution  is 
Gaussian  [44j.  Analysing  such  a  signal  with  continuous  wavelet 
transform  preserves  its  statistical  properties  with  in  the  wavelet 
coefficients.  As  a  result  the  complex  wavelet  coefficients  of  the 
measured  current  i(t)  and  voltage  u(t)  for  a  particular  frequency 
can  be  regarded  as  zero-mean  Gaussian  circular  complex  random 
variables: 

Wi(tJ)  =5R  {Wi(t,f)}+MWi(t,f)} 

Wu(t,f)  =M{Wu(t,f)}+MWu(t,f)},  ™ 

where  Wi(tf)  and  Wu(tf)  are  the  complex  wavelet  coefficients  of 
the  current  i(t)  and  voltage  u(t )  respectively.  Basic  definitions  about 
complex  circular  random  variables  are  given  in  the  Appendix  B. 

Due  to  the  circularity  of  the  complex  random  variables  as  well  as 
the  statistical  properties  of  the  prbs  signal,  the  complex  wavelet 
coefficients  have  the  following  properties: 
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9t{Wu(t,f)},%{Wu(t,f)}  ~  jr(0 yu). 

Using  the  complex  wavelet  coefficients  (5),  the  impedance  at 
particular  frequency /  =  f0  can  be  calculated  as  a  ratio  of  the  com¬ 
plex  wavelet  coefficients  as: 


2.5.  Probability  distribution  function  of  the  impedance  amplitude  jz(  t)  | 

Besides  the  probability  distribution  of  the  real  zr  and  imaginary 
zr  components  of  the  impedance  in  many  cases  the  module  of  the 
complex  impedance  can  be  used  as  a  feature  in  the  process  of 
condition  monitoring.  For  the  circular  Gaussian  complex  random 


z(t) 


Wu(t,f) 

Wi(tJ) 


f=fo 


=  Zr+jZt . 


(6) 


2.3.  Probability  distribution  of  the  complex  impedance  z( t) 

The  ratio  of  two  independent  circular  Gaussian  complex  random 
variables  is  a  known  result  [45],  However,  the  current  i(t)  and 
voltage  u(t)  are  correlated.  This  correlation  can  be  expressed  with 
the  covariance  matrix 

2=[,°"  (7) 

p  a u&i 

where  p  =  pr+jpi  is  complex  correlation  coefficient,  such  as  |p|  <  1. 
As  a  result  the  distribution  of  the  ratio  (6)  reads  [46]: 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

SR (Z)  [O] 


fz(z) 


1  -  \P\2  ( M 2  ,  1  nPrzr-PiZi\ 
rra\o if  ( <rjl  erf  QuCi  J 


(8) 


The  location  of  the  mode  of  fz{z)  depends  on  the  correlation 
coefficient  p.  For  p  =  1/2  +j  1/4  the  shape  of  (8)  is  shown  in  Fig.  2. 


2.4.  Probability  distribution  functions  of  the  real  zr  and  imaginary  Zj 
impedance  components 


The  distribution  of  the  real  zr  and  imaginary  z,  components  of  the 
impedance  (6)  can  be  derived  from  the  distribution  (8).  The  proba¬ 
bility  distribution  functions  of  the  imaginary  component  Zj  reads  [46]: 


fzi(zi\pr,  Pi)  — 


(i  -\p\2y. 


2x°y 


2((1  “  P})ax  +  2 ZiPiOxOy  +  zfafj 


x  3/2 


(9) 


For  the  real  component  zr  the  probability  distribution  function 
has  the  same  shape  and  reads: 


(a)  Distribution  of  5ft(2T) 


(b)  Distribution  of  Q(Z) 


fzr{Zr\Pr,  Pi)  —fzi{Zr\  ~  Pi,  ~Pr)-  (10) 

The  shape  of  the  (9)  and  (10)  are  shown  in  Fig.  3(a)  and  (b). 


Fig.  2.  Contour  plot  of  the  probability  density  function  (8)  for  p  =  1/2  +  jl/4 
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Fig.  3.  Probability  distribution  functions  of  particular  impedance  components  (a)  Distri¬ 
bution  of  $R(Z),  (b)  Distribution  of  §(Z),  (c)  Distribution  function  of  the  module  of  the 
impedance  (12). 
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variables  (5),  their  corresponding  modules  \  Wu(tJ0)  \  and  \  Wi(t.f0)\ 
are  distributed  with  Rayleigh  distribution  as: 

\Wu(t,f0)\  ~  Rayleigh(o-u) 

\Wi(tj0)\  ~  Rayleigh(o-j).  *■  > 

As  the  current  i(t)  and  voltage  u(t )  are  correlated,  the  probability 
distribution  function  of  the  module  of  the  ratio  (6)  reads  [47]: 

(cr?|z|2 +  <r2)  -4p2ala?\z\ 

where  |p|  <  1  is  the  correlation  coefficient,  as  defined  in  (7).  The 
probability  distribution  function  (12)  is  positively  skewed  and 
defined  on  the  non-negative  semi  axis.  A  typical  realisation  of 
the  probability  distribution  function  (12)  is  shown  in  Fig.  3(c). 

2.6.  Parameter  estimation 

The  probability  distribution  functions  (8),  (9)  and  (12)  depend 
on  three  parameters  aUl  a,  and  the  correlation  coefficient  p.  These 
parameters  can  be  easily  estimated  through  the  covariance  ma¬ 
trix  (7)  by  using  the  calculated  wavelet  coefficients  as: 

£{Wu(t,/o)Wu(t,/o)*} 

£{Wi(t,/o)Wi(t,/o)*} 

£{Wi(t,/o)Wu(t,/0)*} 

The  histogram  of  the  impedance  components  and  its  module, 
using  a  single  realisation  of  the  complex  wavelet  coefficients  (5), 
are  shown  in  Fig.  3.  Using  the  estimated  parameters  au,  and  p,  the 
corresponding  probability  distribution  functions  clearly  corre¬ 
spond  to  the  empirically  estimated  histograms. 
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3.1.  Condition  indicator 


Selected  threshold  values  (14)  can  be  additionally  used  for 
mapping  the  complex  impedance  values  into  an  so-called  condition 
indicator  (Cl).  It  is  an  abstract  feature  defined  on  the  interval  [0,oo) 
that  can  be  directly  related  to  the  condition  of  the  system  [35],  For  a 
healthy  system  Cl  should  acquire  low  values,  typically  Cl  =  0.  As  the 
condition  deteriorates,  the  value  of  Cl  should  increase  in  accor¬ 
dance  with  the  severity  of  the  fault(s)  present  in  the  system.  Using 
PFA  value,  the  value  of  Cl  are  typically  scaled  in  such  a  way  that 
values  close  to  one  should  indicate  imminent  failure.  Having  such  a 
feature  value  with  fixed  threshold  determining  imminent  failure 
has  two  major  advantages: 

1.  it  gives  an  unified  and  scale-free  view  of  the  overall  system 
condition  and 

2.  its  trend  can  be  used  for  prognostics  i.e.  estimation  of  the 
remaining  useful  life. 


In  the  context  of  fuel  cells,  the  presence  of  fault  can  either  in¬ 
crease  or  decrease  certain  impedance  components.  Therefore,  for 
fault  free  operation  the  value  of  the  Cl  will  be  in  the  middle  of  the 
region  bounded  by  thresholds  at  certain  PFA,  as  shown  in  Fig.  3. 

The  mapping  into  the  corresponding  Cl  can  be  easily  achieved 
by  using  the  cumulative  distribution  functions  of  either  (9)  or  (12) 
for  the  corresponding  impedance  components.  As  a  result  the 
impedance  components  of  the  fault  free  region  will  acquire  values 
in  the  vicinity  of  1  /2.  Setting  the  maximal  value  for  the  Cl  to  be  the 
chosen  PFA,  the  final  scaling  relation  reads: 


a  =  2  x  (1  -  PFA  r1 


f(z) 


1 

2  ’ 


(15) 


where  F(z)  is  the  cumulative  distribution  function  of  the  corre¬ 
sponding  impedance  component,  as  derived  in  the  Appendix  A. 


4.  Experimental  evaluation 


3.  Fuel  cell  condition  indicator  based  on  the  probability  of 
false  alarm 

The  statistical  properties  of  the  complex  impedance  are 
completely  determined  by  estimating  the  parameters  (13)  using 
the  data  from  the  fault  free  operation.  Deciding  whether  the  fuel 
cell’s  condition  deteriorated  can  be  performed  by  calculating  the 
probability  of  observing  a  particular  impedance  value.  Deter¬ 
mining  the  optimal  values  at  which  the  deterioration  can  be 
regarded  as  significant  is  of  great  importance  for  an  autonomous 
condition  monitoring  system.  The  most  commonly  used  criterion 
for  determining  these  boundary  values  selection  is  low  (PFA) 
[35]. 

Condition  monitoring  system  based  on  low  PFA  will  raise  alarms 
only  when  the  condition  of  the  fuel  cell  is  suboptimal.  Having 
derived  the  probability  distribution  functions  (9)  and  (12),  deter¬ 
mining  optimal  boundary  values  based  on  low  PFA  is  straight  for¬ 
ward  and  can  be  defined  as: 

T  =  1  -  F-1  (PFA)  =  1  -  inf (F(z)  >  PFA},  (14) 

where  F(z)  is  the  cumulative  distribution  function  of  either  (9)  or  ( 12)1 
and  z  is  the  corresponding  impedance  component.  This  approach  was 
used  for  calculating  the  thresholds  values  in  Fig.  3  for  PFA  =  10%. 


1  Corresponding  cumulative  distribution  functions  are  derived  in  the  Appendix  A. 


4.1.  Experimental  setup 

The  test  bed  incorporated  a  PEM  fuel  cell  power  unit,  electronic 
load,  arbitrary  function  generator,  voltage  and  current  measure¬ 
ment  equipment,  and  desktop  computer.  The  measurements  were 
performed  on  a  commercially  available  PEM  fuel  cell  system  HyPM 
HD  8  produced  by  Hydrogenics  Corporation.  The  stack  consists  of 
80  PEM  fuel  cells  each  with  surface  area  of  200  cm2,  providing 
8.5  kW  of  electric  power  in  total.  The  fuel  cell  system  operates  on 
pure  hydrogen  and  ambient  air. 

The  electronic  load  was  used  to  simulate  real  electrical  load  and 
to  perturb  the  fuel  cell  system  with  PRBS  perturbation  signal  in 
galvanostatic  mode.  The  e-load’s  current  was  controlled  by  an 
arbitrary  function  generator,  which  was  generating  the  PRBS 
reference  signal.  Block  diagram  of  the  experimental  setup  is  pre¬ 
sented  in  Fig.  4. 

Voltage  and  current  signals  were  measured  with  an  in-house 
developed  fuel  cell  voltage  monitor  (FCVM)  presented  in  Fig.  5 
[23].  The  top  view  of  the  Fig.  5  shows  the  printed  circuit  board  of 
the  FCVM.  The  board  is  mounted  directly  on  top  of  the  stack  and  the 
galvanic  contacts  to  all  the  fuel  cells  are  implemented  by  spring 
contact  probes,  which  are  visible  in  the  side  view  of  the  Fig.  5.  The 
FCVM  provides  the  means  for  fast  and  accurate  measurements  of 
voltage  and  current  signals  of  any  individual  fuel  cell  inside  the 
stack  (the  contacts  to  all  fuel  cells  of  the  stack  are  clearly  seen  at  the 
bottom  edge  of  Fig.  5).  The  FCVM  utilizes  a  Hall-effect  current 
sensor  to  measure  current,  and  an  analogue  circuitry  to  precisely 
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Fig.  4.  Block  diagram  of  the  experimental  setup. 


measure  changes  in  voltage  of  any  individual  fuel  cell  inside  the 
stack. 

The  FCVM  enables  measurement  of  voltage  and  current  signals 
with  the  resolution  of  80  pV  and  10  mA,  respectively,  at  sampling 
frequency  of  5  kHz.  It  has  a  frequency  range  of  up  to  660  Hz.  The 
FCVM  contains  a  CAN  bus  interface  for  data  transmission  to 
external  device  (e.g.  personal  computer).  Signal  processing  and 
computational  algorithms  were  performed  on  a  PC. 

4.2.  Experimental  procedure 

During  the  experiments,  the  fuel  cell  system  was  kept  at  con¬ 
stant  operating  and  environmental  conditions.  The  fuel  cell  stack 
temperature  and  the  stoichiometry  were  kept  constant  by  the 
built-in  controller  supplied  by  the  fuel  cell  manufacturer  Hydro- 
genics  at  50  °C  and  2.5,  respectively.  At  the  cathode  side,  the 
temperature  and  relative  humidity  of  the  inlet  air  were  measured 
with  Vaisala  HMP235A  transducer  and  controlled  by  external  PID 


Fig.  5.  Top  and  side  view  of  the  FCVM  connected  to  an  8.5  kW  PEM  fuel  cell  system. 


controllers.  The  temperature  of  the  air  flow  was  kept  constant  at 
50  °C  and  the  relative  humidity  changed  in  order  to  introduce 
drying  or  flooding  accordingly.  At  the  anode  side,  the  fuel  cell  was 
fed  with  pure  and  dry  hydrogen  at  constant  temperature  of  20  °C. 
The  DC  current  operating  point  was  set  to  70  A  resulting  in  the 
stack  voltage  of  55  V,  whereas  the  amplitude  of  the  superposed 
PRBS  waveform  was  set  to  2  A.  As  such,  the  peak-to-peak  amplitude 
was  4  A  and  therefore  small  enough  not  to  cause  difficulties  due  to 
non-linearity  of  the  PEM  fuel  cells.  Fig.  6  shows  time  plots  of  cur¬ 
rent  and  voltage  signals  of  a  PEM  fuel  cell  during  measurements. 
Current  was  measured  as  an  absolute  value,  whereas  the  fuel  cell 
voltage  was  measured  as  differential  value. 

During  the  experiment,  an  acquisition  of  measurements 
occurred  every  40  s.  The  experiment  went  through  five  stages.  In 
the  first  stage  of  approx.  11  min,  18  measurements  were  performed 
in  the  fault  free  state  at  which  the  relative  humidity  of  the  inlet  air 
was  kept  at  9%  (which  at  50  °C  corresponds  to  7.5  g  m~3  -  grams  of 
water  per  cubic  metre  of  air).  Then,  air  fed  to  the  fuel  cell  was 
slightly  dried  down  to  relative  humidity  of  5%  (4.1  g  itT3)  for  about 
10  min,  between  the  12  and  the  22  min,  in  order  to  induce  mem¬ 
brane  drying,  and  again  humidified  to  9%.  Afterwards,  between  the 
34  and  the  58  min,  saturated  air  with  relative  humidity  of  100% 
(83  g  m~3),  was  fed  into  the  fuel  cell  system  inducing  flooding  of 
fuel  cells.  Finally,  the  air  humidity  was  decreased  down  to  the  initial 
9%.  Using  such  a  cycle  we  were  able  to  determine  the  response  of 
the  cell  under  different  fault  scenarios  with  various  fault  severity. 
These  stages  of  the  experiment  are  clearly  indicated  in  Fig.  7. 


4.3.  Time  evolution  of  particular  impedance  components 

Using  the  complex  wavelet  coefficients  (4),  calculated  from  the 
initial  fault  free  operation  signals,  the  parameters  of  probability 
distribution  functions  (8)  and  (12)  were  estimated  employing  (13). 
With  the  estimated  parameters,  the  threshold  values  were  calcu¬ 
lated  for  four  different  PFA  =  {25%,  20%,  15%,  10%}.  The  resulting 
thresholds  for  each  impedance  component  are  shown  in  Fig.  7.  The 
threshold  values  clearly  define  the  fault  free  region. 

The  first  departure  from  the  fault  free  condition  was  between 
the  12th  and  the  22nd  minute,  when  the  fuel  cell  was  fed  with 
somewhat  drier  air  (5%  relative  humidity).  In  that  particular  region 
the  real  component  of  the  impedance  increased  and  surpasses  the 
threshold  of  70%,  as  shown  in  Fig.  7(a)  and  (b).  Similar  observation 
can  be  made  also  for  the  module  of  the  impedance,  as  shown  in 
Fig.  7(c).  Despite  the  minute  changes  in  the  relative  humidity  of  the 
inflow  air,  there  are  clear  change  in  the  value  of  the  impedance, 
with  30%  probability  of  false  alarm.  Without  the  threshold  values 
such  an  quantification  would  have  been  impossible. 
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Fig.  6.  Measured  current  and  voltage  signals  during  the  experiment. 
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Fig.  7.  Trend  of  impedance  components  throughout  the  experiment  (a)  Values  of  3 J(Z), 
(b)  Values  of  3(Z),  (c)  Values  of  |Z|. 

Unlike  the  minor  changes  caused  by  the  mild  decrease  in  the 
humidity,  the  flooding  with  100%  saturated  air  inflow  shows  sig¬ 
nificant  changes  in  al  three  impedance  components.  The  values  of 
Tt(z)  and  |z|,  between  the  34th  and  the  58th  minute,  indicate  that 
the  condition  of  the  fuel  cell  departed  from  the  fault  free  area,  with 
less  then  10%  probability  of  false  alarm,  as  shown  in  Fig.  7(a)  and  (c). 
The  change  in  the  imaginary  component  3(z)  is  somewhat  smaller, 
the  probability  of  false  alarm  is  less  then  30%. 


Finally,  when  the  relative  humidity  of  the  air  inflow  was 
returned  to  normal  the  impedance  components  returned  to  the 
initial  values  within  the  fault  free  region. 

4.4.  Time  evolution  of  the  Cl 

The  time  evolution  of  the  Cl  (15)  of  the  impedance  module  |z|  is 
shown  in  Fig.  8.  The  values  of  the  Cl  were  calculated  using  the 
mapping  relation  (15)  with  PFA  =  90%. 

By  analysing  this  time  evolution  one  can  notice  that  during  the 
fault  free  operation  Cl  acquires  low  values.  As  the  first  fault  was 
induced,  the  values  of  Cl  increased  up  to  0.5,  which  indicates  that 
the  condition  of  the  fuel  cell  deteriorated.  The  acquired  value  of  0.5 
clearly  shows  departure  from  fault  free  operation  while  in  the  same 
time  is  inline  with  the  low  severity  of  the  induced  drying,  relative 
humidity  change  of  8%.  Unlike  the  low  severity  of  the  drying  fault, 
during  the  flooding  phase  of  the  experiment,  the  value  of  Cl  sur¬ 
passes  the  limit  value  of  1.  The  overwhelming  flooding  of  the  sys¬ 
tem  is  directly  reflected  in  the  high  value  of  the  Cl.  Finally,  at  the 
end  of  the  experiment,  when  the  air  condition  was  back  to  nominal, 
the  values  of  Cl  decreased  indicating  almost  fault  free  operation. 

The  results  confirm  the  behaviour  of  Cl.  It  is  sufficiently  sensitive 
on  deviation  in  the  condition  of  the  water  management.  In  the 
same  time  its  value  is  directly  related  to  the  severity  of  the  devia¬ 
tion  with  respect  to  the  selected  PFA. 

5.  Conclusion 

Observing  the  impedance  of  a  fuel  cell  in  the  context  of  random 
variables,  this  paper  specifies  probability  distribution  functions  of 
the  real  and  imaginary  components  as  well  as  its  module.  Using  the 
data  solely  from  the  fault  free  operation,  one  can  estimate  the  free 
parameters  of  these  distribution  functions.  As  a  result,  threshold 
values  defining  the  fault  free  region  can  be  calculated  by  specifying 
a  specific  probability  of  false  alarm. 

Determining  the  threshold  values  using  the  probability  of  false 
alarm  has  significant  practical  merit.  For  each  of  the  measured 
impedance  components  one  can  calculate  the  probability  of  false 
alarm  and  use  this  value  as  an  indicator  of  confidence  that  the 
current  condition  of  the  fuel  cell  has  sufficiently  departed  from  the 
fault  free  region.  Furthermore,  the  proposed  condition  indicator 
scaled  with  the  chosen  PFA  is  directly  related  to  the  fault  severity 
and  can  be  used  for  both  fault  detection  as  well  as  estimation  of  the 
remaining  useful  life.  The  evaluation  of  the  proposed  method 
confirms  that  the  calculated  condition  indicator  based  on  the  PFA 
corresponds  to  the  severity  of  the  induced  faults. 


Fig.  8.  Trend  of  Cl  of  |Z|  throughout  the  experiment. 
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Appendix  A.  Cumulative  distribution  functions 


^[y2]  (B.3) 

Therefore  the  random  variables  X  and  Y  are  independent  with 
same  variance. 

Appendix  B.l  Circularity  of  spectral  components 

A  stationary  random  signal  can  be  represented  using  Cramer 
representation  [49]: 


£[X2] 

E[XY] 


Appendix  A.i.  Cumulative  distribution  function  of  Zj 

According  to  46]  integrating  out  zr  from  (8)  results  into  the 
following  cumulative  distribution  function: 
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Solution  of  (A.I )  for  arbitrary  y  =  PFA  reads: 
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Appendix  A.2.  Cumulative  distribution  function  of  \z\ 

The  cumulative  distribution  function  that  corresponds  to  (12) 
reads  [47]: 


F|z|(z): 


2\jAz2 0  -p2)afai+  (°l 


(A.3) 


where  z  >  0.  Solution  of  (A.3)  for  y  =  PFA  reads: 


CO 

x(t)  =  J  e,‘jtdZ(w), 

— oo 


(B.4) 


where  dZ(w)  are  complex  random  variables  directly  related  to  the 
signals  power  spectrum  with  the  following  properties: 

E[dZ(oj )]  =  0,  £[dZ(w1)dZ(w2)]  =  0,  for  oj^cj2.  (B.5) 

Therefore  for  stationary  signals  the  complex  random  variables 
dZ(w)  are  circular  [48], 
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A  Gaussian  complex  random  variable  can  analysed  through  its 
real  and  imaginary  components 
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